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EXECUTIVE  SUMMARY 


Of  the  earth  filters  that  determine  seismic  waveforms, 
the  Q  filter  is  probably  the  least  well  constrained.  Though 
anelastic  attenuation  reduces  signal  amplitude  by  several 
orders  of  magnitude  at  5  Hz,  the  form  of  the  Q  filter  is 
still  in  question.  In  this  study,  a  form  is  adopted  for 
Q ( f )  which  meets  the  constraints  of  both  microphysics  and 
continuum  mechanics.  This  is  turn  leads  to  the  definition 
of  t*(f)  in  terms  of  x^,  a  path  averaged  absorption  band 
parameter.  With  the  completion  of  the  forward  problem,  an 
inverse  theory  is  completed  with  the  capability  to  model 
several  forms  of  the  Q  filter  and  to  estimate  parameters  of 
each  model. 

The  inverse  process  is  applied  to  two  spectral  ratio 
data  sets,  PP/P  and  receiver  function  ratios.  The  PP/P  data 
were  strictly  long  period,  so  only  the  difference,  6t*,  in 
low  frequency  t*  values  could  be  estimated.  In  practice, 
the  Q  filter  signal  is  masked  by  noise  on  the  spectrum  due 
to  the  crustal  transfer  function  at  the  PP  bounce  point. 
The  inverse  method  is  demonstrated,  but  adequate  estimates 
of  attenuation  properties  will  require  broadband  seismic 
data . 

The  receiver  function  spectral  ratios  are  short  period 

and  thus  provide  much  better  resolution  on  Q.  Consistent 

values  of  fit*  are  found  for  several  paths  about  the  Eastern 

Kazakh  and  NTS  nuclear  test  sites.  In  tests  of  model  def- 
.  .  * 

mition,  the  frequency  dependent  t  models  are  shown  to  be 
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significantly  better  at  fitting  high  frequency  data  than  the 

•k 

constant  Q  (or  constant  t  )  assumption.  Bandwidth  and 
instrumental  gain  problems  again  prevent  any  absolute  con¬ 
clusions  . 
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1.  INTRODUCTION 

The  anelastic  attenuation  of  a  seismic  wave  can  be 
*  - 1 

written  in  terms  of  t  ,  the  weighted  average  of  Q  over  the 

k 

seismic  raypath.  But,  because  t.  is  a  path  property,  esti- 

* 

mation  of  t  requires  separation  of  anelastic  effects  from 
source  and  geometric  properties.  Spectral  ratios  have  been 
very  successful  in  such  separation,  using  the  ratio  of 
seismogram  at  one  receiver  to  that  at  another  or  the  ratio 

■k 

of  two  phases  at  the  same  receiver  to  estimate  6t  ,  the 

k  k 

difference  in  t  for  the  two  paths.  Typically  t  is  assumed 
to  be  independent  of  frequency,  and  (St  is  estimated  from  a 
linear  fit  to  the  graph  of  log  spectral  ratio-vs- frequency . 
Strong  evidence  exists,  however,  that  Q  and  therefore 

k 

t  must  be  frequency  dependent  (Lundquist  and  Cormier, 
1980;  Lundquist  et  al,  1980),  and  that  the  fit  of  a  straight 
line  is  therefore  not  appropriate.  This  report  will  discuss 

k 

the  choice  of  a  model  for  frequency  dependent  t  and  will 

* 

develop  the  framework  for  inverting  for  t  without  the 
restrictive  assumption  of  frequency  independence.  The 
inversion  technique  will  then  be  tested  against  two  dif¬ 
ferent  types  of  spectral  ratios. 


I 


I 


SG1-R-81-052 


‘I 


t 


f 


♦ 


♦ 


♦ 


♦ 


I 


1 1 .  THEORY 

The  development  of  an  inverse  theory  for  frequency 

* 

dependent  t  requires  three  steps.  First,  the  foreward 
problem  must  be  solved  for  seismic  raypaths  in  general. 
Second,  the  specific  spectral  ratios  under  consideration 
must  be  examined  and  understood  in  detail.  Finally,  the 
inverse  theory  must  be  cast  in  terms  of  the  relative  be¬ 
havior  of  parameters  in  the  foreward  problem.  Each  of  these 
steps  will  be  discussed  in  a  separate  section. 

2,1,  JSLlfJ. 

Q  models  may  be  constrained  a  priori  to  certain  classes 
of  frequency  dependence  and  certain  forms  of  depth  depend¬ 
ence  (Lundquist,  1979).  An  examination  of  anelasticity  from 
both  the  continuum  and  microphysics  points  of  view  shows 
that  all  of  the  important  physical  dissipation  mechanisms 
have  two  characteristics  in  common  (Jackson  and  Anderson, 
1970).  They  are  all  thermally  activated  processes,  and  they 
are  all  relaxation  processes.  Thus  a  general  form  may  be 
written  for  the  frequency  dependency  of  a  single  mechanism 
without  adopting  any  particular  mechanism  for  attenuation  in 
the  earth.  The  of  this  standard  anelastic  solid  is  a 

perked  function  of  frequency  about  w  =  1/t ,  where  t  is  the 
relaxation  time. 

A  distribution  of  such  mechanisms  may  be  assumed  to 
exist  in  a  complex  solid  such  as  the  earth,  giving  at¬ 
tenuation  over  a  broad  range  of  frequencies.  Liu  et  al 
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(1976)  adopted  the  distribution  in  terms  of  relaxation  time 
given  by  Gurevich  (1964): 

1/x,  x1  <  r  <  t2  sec/rad  (1) 

0  ,  otherwise 

This  defines  an  absorption  band,  as  shown  in  Figure  1,  where 
and  t2  are  the  half  amplitude  points  at  low  and  high 
frequencies,  respectively.  An  absorption  band  rheology  is 
clearly  only  a  first  approximation  to  the  true  frequency 
dependence  of  Q,  but  it  is  the  simplest  Q  (f)  model  compat¬ 
ible  with  the  microphysics  of  dissipation  mechanisms. 

The  significance  of  activated  process  theory  for  Q(f,r) 
is  that  the  theoretical  depth  dependence  of  relaxation  time 
is  an  exponential  function  of  pressure  and  inverse  tem¬ 
perature.  Essentially,  t  increases  rapidly  with  depth 
through  the  range  of  the  asthenosphere ,  where  increasing 
temperature  dominates.  Then  x  decreases  slowly  as  pressure 
effects  dominate.  The  result  is  that  the  absorption  band 
for  a  given  mechanism  must  shift  to  high  frequency  over  the 
asthenosphere  and  then  to  lower  frequency  with  increasing 
depth.  If  the  same  mechanism  is  not  active  at  all  depths, 
then  the  picture  becomes  much  more  complex. 

The  conclusion,  then,  is  that  although  the  assumption 
thrt  Q  1  may  be  modeled  by  an  absorption  band  is  valid  on  a 
pointwise  basis,  it  is  only  approximately  correct  on  a 
pathwise  basis.  The  assumption  of  constant  i2  is  the  same  as 
saying  that  Q-1  (3Hz),  say,  is  the  same  fraction  of  Qm_1  at 


D(t)  = 


SGI-R-81-052 


ABSORPTION  BAND 


1.  Normalized  absorption  band  versus  frequenc 
t \  and  T2  are  relaxation  times  at  the  low- 
and  high-  frequency  half-amplitude  points, 
respectively. 
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all  depths.  Because  a  seismogram  is  a  path  average,  then 
any  estimate  of  t 2  measured  from  seismogram  spectral  ratios 
is  a  weighted  average  of  t2  over  the  path.  All  values  of  t0 

reported  in  this  paper  will  have  this  restricted  definition. 

* 

t  is  a  path  property  which  depends  upon  a  weighted 
average  of  Q-1  over  the  path.  The  adoption  of  an  absorption 
band  model  for  a  uniform  medium  leads  to  the  definition. 


*  * 
t  (u>)  =  t 

m 


_2 

n 


tan 


-1 


iu(t  -t  ) 

1  2 

1+tU^T  T 
1  2 


which  reduces  to 


t*(u>) 


(2) 


(3) 


.  *  * 
under  the  assumption  that  t  is  the  t  at  the 

"k 

maximum  of  the  absorption  band,  corresponding  to  t  values 
reported  in  the  literature  for  low  frequencies.  The  path  Q 

: k 

filter  is  given  in  terms  of  t  by 


ln|FQ|  =  -»f  t*m 
where  f  is  frequency. 


(4) 


The  spectral  ratio  of  two  seismograms  can  be  defined  in 
terms  of  equation  (4)  by 
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» 

t  1 


log  ZA  =  -2fB 
S 

o 


t* 

mO 


(5) 


where  C  is  a  relative  gain  factor.  The  spectral  ratio  is 

.  .  *  * 
defined  m  terms  of  five  parameters:  t  A  t™,  t  ,  ,  , 

mi),  20  mk  2k 

and  C.  Since  the  ratio  is  determined  at  many  frequencies, 
the  inverse  problem  becomes  a  least  squares  process. 
Equation  (5)  will  be  the  basis  for  the  inverse  theory  pre¬ 
sented  here. 

Figures  2  and  3  show  several  examples  of  theoretical  Q 
filters,  corresponding  to  equation  (4),  and  synthetic  spec¬ 
tral  ratios  as  defined  by  equation  (5).  A  feature  of  Q 

* 

filters  defined  by  a  frequency  dependent  t  is  that  the 
amplitude  spectrum  on  a  log-log  plot  has  two  inflections, 

■k 

one  determined  by  t  and  one  by  t For  t2,  greater  than 
zero,  the  filter  approaches  a  constant  spectral  level  at 
high  frequencies,  in  direct  contrast  to  the  frequency  in¬ 
dependent  filters  which  fall  exponentially  at  high  fre¬ 
quency.  This  difference  can  amount  to  orders  of  magnitude 
difference  in  predicted  spectral  amplitudes,  and  will  be 
significant  in  determining  whether  is  a  required  para¬ 
meter  in  the  Q  filter. 


2.2  Spectral  Ratios 

Two  specific  spectral  ratio  types  will  be  discussed  and 
tested  in  this  report.  The  ratio  of  PP  waves  to  the  direct 

•k 

P  arrival  produces  a  6t  which  is  a  direct  estimate  of  the 

* 

t  at  the  bounce  point.  The  high  quality  spectral  ratios 
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» 
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defined  during  estimation  of  relative  receiver  functions 

* 

(Lundquist  et  al,  1980)  produce  6t  which  represents  the 
difference  in  anelastic  attenuation  between  the  two  seismic 
raypaths.  These  ratios  use  the  approximately  symmetric 
nuclear  explosion  sources,  and  several  estimates  of  the 
ratio  are  averaged  to  minimize  contamination  by  source 
asymmetries . 

A  P  wave  seismogram  trace  may  be  defined  by 

Sp  (t)  =  F(t)  *  Cs  *  Gp  *  Qp  *  CR  *  I  (6) 

where  F  is  the  source  function;  Cg  and  Cr  are  crustal  trans¬ 
fer  functions  at  source  and  receiver,  respectively;  G  is 
the  geometric  attentuatiori;  I  is  the  instrument,  and  Qp  is 
the  desired  Q  filter.  The  corresponding  PP  wave  is 

Spt(t)  =  H  |F(t)|  *  Cs  *  Cpp  .  Qpp  *  CB  *  CB  *  Dpp  *  CR  *  I 

where  H  implies  Hilbert  transform;  CR  is  the  crustal  trans¬ 
fer  function  at  the  bounce  point,  and  D  is  the  surface 
reflection  coefficient.  The  Hilbert  transform  affects  only 
the  phase  and  is  necessary  for  synthesis  of  PP  from  P,  but 
it  is  not  a  necessary  correction  for  spectral  ratios. 

The  spectral  ratio  then  is 

R(u) )  =  Spp  ( u ) )  =  Gpp  •  Dpp  •  Qpp  •  CB  •  CB 

S  (ui)  G  Q 

P  P  P 
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Here  the  factors  have  been  grouped  to  isolate  the  frequency 
independent  quantities  G  and  D  from  those  with  frequency 
dependence.  The  assumption  that  the  source  factors  divide 
out  is  a  statement  that  take-off  angles  do  not  vary  signifi¬ 
cantly  between  P  and  PP  rays.  For  shallow  focus  and  dis¬ 
tances  of  80°  and  40°,  this  implies  an  assumption  that  the 
source  and  source  crustal  response  do  not  change  over  about 
10°  change  in  emergence  angle.  Violation  of  this  assumption 
will  be  seen  as  a  noise  factor  (or  variance)  in  the  spectral 
ratio. 

In  terms  of  log  amplitude  spectra, 

J  R  (  f )  |  =  -nfB  ^  t*pp"t*p  ^  +  log  ^-PP.-PP^  +  2B  log  CB(f)  (9) 


Equation  (9)  defines  the  basic  data  set  and  foreward  solu¬ 
tion  for  P/PP  spectral  ratios.  The  correspondence  with 
equation  (5)  is  clear,  except  that  a  noise  term  has  now  been 
added  in  terms  of  the  frequency  dependent  bounce  point 
crustal  response.  In  practice,  the  spectral  ratio  will  be 
smoothed  by  a  running  average  to  minimize  the  effects  of  Cg. 
That  is,  the  Q  filter  is  expected  to  cause  the  broad  fre¬ 
quency  dependence  in  the  spectral  ratio,  while  the  crustal 
response  is  expected  to  cause  a  variance  about  the  desired 
signal . 

A  A 

The  fit  determined  by  (9)  may  be  interpreted  as  the  t 
associated  with  the  crust  and  upper  mantle  at  the  bounce 
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point,  the  only  significant  difference  between  the  two 
travel  paths.  Note  that 


(A) 


2  t*  (A/2) 


But  t* 

P 


has  been  observed  to  be  nearly  constant  with  dis¬ 


tance,  so 


*  *  * 
tpp“tP-  tP 

which  may  be  interpreted  directly  as  a  property  of  the 
bounce  -ioint  Indeed,  this  method  may  provide  the  only 
correct  definition  of  t  as  a  function  of  geographic  loca¬ 
tion  . 

U  is  interesting  to  observe  that  each  of  the  terms  in 
(9)  may  eventually  be  recovered.  The  reflection  coefficient 
and  geometric  spreading  factors  can  be  estimated  from  stand¬ 
ard  models,  and  any  error  will  be  recovered  from  the  gain 
factor,  C,  in  (5).  This  is  true  in  the  sense  that  the 
reflection  coefficient  is  actually  the  zero  frequency  re¬ 
sponse  of  the  crust.  Since  (uj)  may  be  assumed  to  have  no 
broad  frequency  trends,  then  the  ratio  corrected  for  the 
estimated  6t  is  a  direct  estimate  of  the  amplitude  spectrum 
of  CB. 

As  a  final  note  on  PP/P  spectral  ratios,  an  independent 
estimate  of  the  quality  of  the  St  values  may  be  obtained 
from  time  domain  synthesis  of  the  PP  waveform  from  the  P 
waveform. 
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spp  <“»  = 


H 


{sp(“)}  •  s  («t‘) 


_££ 


D_ 


PP 


(10) 


The  error  between  the  estimate,  Sp^  and  the  observed  PP 
seismogram  will  result  primarily  from  inability  to  model  C^, 
but  the  comparison  of  waveforms  will  be  a  verification  of 
the  assumptions  made  in  this  analysis.  Poor  time  domain 

4 

fits  can  also  be  an  indication  of  the  presence  of  uncor¬ 
related  high  frequency  noise,  which  is  treated  as  signal  by 
spectral  ratio  technique. 

Receiver  function  spectral  ratios  differ  from  PP/P 
ratios  in  that  seismograms  from  two  separate  stations  are 
used  in  the  ratio.  Thus,  no  factors  may  be  assumed  a  priori 
to  divide  out. 


R(u>) 


=  Gi  .  Q  i 
S  (w)  G  Q 
o  o  o 


(11) 


Here  the  factors  have  been  grouped  according  to  importance . 
The  ratio  of  geometric  spreading  factors  is  a  simple  gain 
factor  with  no  frequency  dependence.  The  factors  in  the 
last  bracket  will  divide  out  if  the  source  is  symmetric  (an 
explosion),  the  source  crustal  layering  is  flat,  and  the 
instrument  responses  are  known.  Unfortunately,  though  the 
shape  of  the  instrument  response  curves  are  known,  gain 
factors  often  are  not.  Ideally,  both  G  and  I  are  known;  so 
that  C  in  (5)  is  zero.  In  practice,  azimuthal  and  take-off 
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» 


angle  asymmetries  exist  in  S  and  C  ,  leaving  uncancelled 
frequency  dependent  noise,  N(w),  in  the  spectral  ratio.  But 
the  asymmetries-  are  dependent  upon  source  location  and 
source  depth,  so  averaging  over  several  sources  confined  to 
a  limited  region  (an  explosion  test  site)  reduces  the  noise 
factor.  This  process  produced  clean  relative  receiver 
functions  for  stations  about  the  Eastern  Kazakh  test  site  in 
Russia  (Lundquist  et  al,  1980c),  but  was  less  successful  for 
the  more  complex  NTS  in  Nevada  (Lundquist  et  al,  1981). 

With  these  observations  in  mind,  (11)  can  be  reduced  to 


log  R(f) 


-n  fB( t*-t* )  +  B  log  C  +  B  log 


'Ri 


'Ro 


+  N(f) 


(12) 


t 


* 


where  C  includes  geometric  and  instrumental  gain  factors. 
For  Q  inversion,  the  ratio  of  receiver  crustal  responses 
becomes  the  most  important  noise  term,  dominating  the  re¬ 
sidual  noise  due  to  source  effects.  The  variance  due  to 
these  factors  can  to  some  extent  be  reduced  by  smoothing  the 
spectral  ratio,  at  which  time  equation  (12)  defines  the  data 
set  and  foreward  problem  for  the  receiver  function  spectral 
ratios . 


2.3  The  Inverse  Problem 

The  inverse  problem  for  a  frequency  dependent  t*  can  be 
written  as  an  error  minimization  problem  with  five  unknowns 
in  the  presence  of  noise.  This  is  true,  under  the 


r. 


L 
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assumptions  stated  here,  for  a  wide  range  of  spectral  ratio 
types.  Thus  an  inverse  program  will  be  a  useful  tool  for 
the  analysis  of  path  attenuation  in  a  variety  of  circum¬ 
stances  . 

Let  the  solution  the  the  foreward  problem  define  a 
vector 

g(u>)  =  c  h  (f,  tQ*,  t 0 ,  tR*,  t k )  (13) 

where  c  is  the  overall  gain  factor,  and  h  implies  function. 
Using  the  notation  g^  =  g  ( f ^ ) ,  the  partial  derivative 
matrix  may  be  defined  as: 


=  ^  =  c'  ^  ’i'  t>2'  >2 

Then  an  increment  of  change  in  the  model  parameters  produces 
an  increment  of  change  in  g  according  to 


The  observed  spectral  ratio,  r^,  may  be  compared  to  the 
foreward  problem  solution  to  define  an  error  vector 


di s  u  -  «i 


The  object  of  the  inversion  is  to  find  a  set  of  m^  such  that 
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the  RMS  quantity 

IK  '  Aij  6mjll 

is  minimized.  The  solution  for  6m^  gives 

6m-  =  (ATA)_1  AT  d  (14) 

where  the  ~  implies  estimate. 

In  the  presence  of  noise 

6mj  =  (ATA  +  o2I  )_1  AT  d  (15) 

where  o  is  the  noise  variance.  If  weights  are  added, 

6rhj  =  (AT  WA  +  o 2 1  )  ~ 1  AT  Wd  (16) 

The  partial  derivatives  are  given  by 


=  a 


3mj  6mj 
9gi 


at. 


-2fB 


*  =  -2fB 


t*  tan'1/  — 


tan 


-1 


UJ  T 

v  i  k 


U)T 


O  J 


-  t 


*tan 


-1 


w  i 

i  o 


+C 


3g 


3  T 


=  -2fB 


* 

t  uj  • 
o  1 


,.22 

1+UJ  T 

o 


dJh 

at* 


-2fB  [tan-1  1 

L  ujt  k  J 
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-2fB 


* 

t  10 


i  k 


3g 

_ i  =  1 

3c 


Figure  4  shows  an  inversion  run  on  noise-free  synthetic 
data.  The  computer  caption  underneath  the  figure  provides 
the  critical  information  on  this  run,  including  final  values 
for  the  five  variables,  the  starting  values,  weights,  number 
of  iterations  and  the  error.  The  bandlimited  data  set  used 
is  representative  of  the  widest  frequency  range  likely  to  be 
realized  in  short-period  spectral  ratios  of  digital  data. 
Bandlimiting  severly  reduces  the  resolution  in  the  inverse, 
and  is  partially  responsible  for  the  error  in  the  final 
values  of  model  parameters. 

Another  source  of  error  in  these  inversions  is  trade- 

*  * 
off,  both  between  t  and  and  between  the  two  t  '  5  and  C. 

Thus  the  starting  values  and  judicious  use  of  weights  is 

critical  to  obtaining  correct  results.  In  the  work  to  be 

presented  in  the  next  section,  multiple  tests,  reduction  of 

variables  and  solution  for  differences  rather  than  absolute 

values  will  all  be  used  to  reduce  the  effects  of  tradeoff. 
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1  I_I_ DATA  AND  RESULTS 

3.1  PP.  P  Spectral  Ratios 

The  spectral  ratio  of  PP/F  waves  may  be  interpreted  in 
terms  of  anelastic  properties  in  the  crust  and  upper  mantle 
of  the  bounce  point.  Two  bounce  points  were  chosen  for 
study  here,  one  at  the  Soviet  nuclear  test  site  in  eastern 
Kazakhistan,  and  one  in  the  east  Pacific  near  the  Hawaiian 
Islands.  The  first  represents  conditions  of  an  old  thick 
continental  crust,  while  the  second  represents  a  young  thin 
oceanic  crust.  A  third  bounce  point  at  NTS  representing  a 
tectonic  crust  will  be  examined  in  a  future  report. 

For  each  bounce  point,  source  receiver  pairs  were 
sought  which  had  distances  in  the  range  60  to  85° .  For  the 
Soviet  bounce  point,  sources  in  the  Kuriles  and  Japan  were 
recorded  in  the  Mediteranean  region  and  sources  in  southern 
Asia  were  recorded  in  Scandinavia.  This  provided  a  wide 
azimuthal  coverage,  as  is  shown  in  Figure  5,  around  a  region 
of  about  5°  radius  (550km)  from  the  Russian  nuclear  test 
sites.  The  oceanic  crust  was  studied  using  Tonga-Fiji 
earthquakes  recorded  at  stations  in  the  western  U.S.  These 
bounce  points  were  constrained  to  be  further  than  2°  from 
the  Hawaiian  Islands,  so  that  the  reflection  would  occur  in 
simple  crust.  Table  1  lists  the  source/receiver  pairs  used 
in  this  study. 

WWSSN  data  were  used  for  this  part  of  the  study  because 
of  the  relative  rarity  of  paths  to  digital  stations.  Though 
short  period  records  were  sought,  either  the  PP  was  below 
noise  level  or  the  P  was  off  scale.  Thus  this  analysis  is 
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TABLE  1 

Data  set 

for  PP/P  t*  inversion 

KAZAKH  BOUNCE  POINT 

REGION 

DATE 

DEPTH  A 

STATION 

Near  E.  Coast  Honshu, 

Japan 

3-16-65 

34.00  85.41° 

HLW 

Ryukyu  Islands 

9-21-65 

199.00  82.9° 

ATU 

So.  Honshu,  Japan 

8-13-67 

357.00  83.41° 

ATU 

Ryukyu  Islands 

5-14-68 

162.00  83.16° 

ATU 

Shikoku,  Japan 

8-5-68 

41.00  82.81° 

ATU 

Sea  ol  Japan 

3-31-69 

417.00  81.43° 

HLW 

E.  Russia-N.E.  China 

Border 

4-10-69 

555.00  73.39° 

HLW 

77.18° 

JER 

Northeastern  China 

7-18-69 

33.00  66.48° 

1ST 

Burma -India  Border 

10-17-69 

134.00  60.88° 

KEV 

Sakha  1 i n  Island 

12-18-69 

344.00  80.13° 

EIL 

Eastern  Sea  of  Japan 

6-10-71 

226.00  82.43° 

HLW 

E.  Russia-N.E.  China 
Borde  r 

9-10-73 

532.00  74.75° 

EIL 

Hokkaido,  Japan 

11-8-74 

130.00  83.75° 

HLW 

Northeastern  China 

2-4-75 

33.00  67.14° 

1ST 

PACIFIC  OCEAN  BOUNCE  POINT 

REGION  DATE 

DEPTH 

A 

STATION 

Fiji  1-26-72 

668.00 

77.10° 

BKS 

Fiji  5-9-72 

568.00 

76.74° 

GSC 

Fiji  12-19-73 

246.00 

79.51° 

COR 

77.83° 

GSC 

Fiji  7-6-77 

594.00 

78.16° 

BKS 
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limited  to  long  period  WWSSN  seismograms.  Each  seismogram 
was  digitized  four  times,  and  each  digital  copy  was  compared 
to  the  original  to  verify  the  accuracy.  The  four  traces 
were  then  averaged,  corrected  for  paper  magnification,  skew 
of  the  recording  system  and  instrument  gain  factor. 

The  resulting  spectral  ratios  had  a  useful  bandwidth  of 
one  decade,  from  .05  to  .5  Hz.  Because  the  effects  of 
anelasticity  is  small  at  low  frequencies,  and  because  x2  has 
been  found  to  be  a  high  frequency  parameter  (Lundquist  and 
Cormier,  1979),  the  parameters  in  the  PP/P  inversions  were 

•k  k 

reduced  to  t  ^  11  pp  ancl  C'  by  setting  x2  to  zero  and  ap¬ 

plying  weights  to  keep  it  small .  Without  frequency  de- 

* 

pendence,  there  is  total  tradeoff  between  the  t  values,  so 

the  results  of  this  experiment  can  only  be  interpreted  in 
* 

terms  of  6t  m,  the  low  frequency  difference  in  average  path 
attenuation. 

An  example  of  a  PP/P  spectral  ratio  inversion  is  shown 
in  Figure  6  for  a  Kurile  source  recorded  at  station  EIL  in 
Israel.  Both  time  and  frequency  domain  results  are  shown, 
and  both  waveform  fitting  and  spectral  fitting  are  com- 

k 

patible  with  6t  =  0.8sec.  The  amplitude  of  the  variance  on 
the  spectrum  is  small,  implying  a  simple  crust  at  the  bounce 
point.  The  time  domain  results  confirm  this  conclusion 
showing  that  the  PP  waveform  is  the  theoretical  Hilbert 
transform  of  the  P  wave. 

Unfortunately,  the  success  of  Figure  6  is  not  typical. 
The  range  of  scatter  is  given  as  a  histogram  in  Figure  7, 


I 
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Figure  6.  PP/P  spectral  ratio  inversion  for  9/10/73  recorded 
at  station  EIL  with  6t*=0.8sec.  For  the  waveform 
comparison,  the  observed  PP  wave  is  shown  by  the 
solid  line,  the  synthetic  PP  is  shown  by  the  dashed 
line. 
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and  tour  additional  examples  are  shown  in  Figures  8-11  with 
* 

a  range  in  fit  of  .27  to  1.4  sec.  Examination  of  the  spec¬ 
tra  shows  that  each  one  has  a  much  larger  variance  than  that 
shown  in  Figure  6.  The  conclusion  must  be  that  the  noise 
added  to  the  spectrum  at  the  crustal  bounce  point  is  com¬ 
monly  so  great  as  to  mask  the  anelastic  attenuation  seen  by 
these  low  frequencies. 

Two  simple  tests  are  shown  in  Figure  12  to  demonstrate 
the  level  of  signal  sought  by  the  PP/P  inversion.  The  plot 

k 

on  the  light,  shows  theoretical  spectra  at  6 1  values  of  .4 

to  1.2  sec.  For  data  bandlimited  at  .6  Hz.,  the  maximum 

* 

amplitude  variation  for  a  difference  in  6t  of  0.8  is  a 
factor  of  four,  and  most  of  the  significant  amplitude  '’illa¬ 
tion  occurs  in  the  range  .3  to  .  5Hz .  Thus  noise  at  these 

frequencies  can  easily  dominate  the  effects  of  ane lastici ty . 

* 

The  time  domain  test  confirm  that  a  change  of  .8  sec  in  fit 
has  virtually  no  effect  on  the  results.  That  is,  there  is  no 

k 

resolution  in  fit  using  long  period  time  domain  synthesis. 

Similar  results  are  observed  for  the  Pacific  Ocean 

bounce  point,  as  is  shown  in  Figures  13  and  14.  These  tests 

found  that  rather  than  being  simpler,  the  oceanic  bounce 

point  produces  a  more  complex  waveform.  The  example  in 

Figure  13  is  one  of  the  few  cases  where  the  PP  is  a  simple 

Hilbert  transform  of  the  P  wave.  In  addition,  a  number  of 

PP  waves  were  observed  to  have  greater  amplitude  than  the 

* 

coi  respond! ng  P  wave,  giving  a  negative  fit  .  Such  a  con- 
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Observed  P  wave 
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Figure  8.  PP/P  spectral  ratio  inversion  for  2/4/75  recorded 
at  station  1ST  with  6t*=0.27sec.  For  the  waveform 
comparison,  the  observed  PP  wave  is  shown  by  the 
solid  line,  the  synthetic  PP  is  shown  by  the  dashed 
line. 
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Figure  9.  PP/P  spectral  ratio  inversion  for  8/13/67  recorded 
at  station  ATU  with  St*=0.55sec.  For  the  waveform 
comparison,  the  observed  PP  wave  is  shown  by  the 
solid  line,  the  synthetic  PP  is  shown  by  the  dashed 
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Figure  10.  PP/P  spectral  ratio  inversion  for  4/10/69 
recorded  at  station  JER  with  St*=0.99sec. 
For  the  waveform  comparison,  the  observed 
PP  wave  is  shown  by  the  solid  line,  the  syn 
thetic  PP  is  shown  by  the  dashed  line. 
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Figure  12.  Test  results  are  shown  for  synthetic  seismograms  and 

Q  filters  as  a  function  of  fit*.  The  observed  seismogram 
(dashed  line)  is  compared  with  five  different  synthetic 
seismograms  (solid  line)  computed  with  different  fit* 
values.  Three  different  Q  filters  are  compared  to  show 
their  effect. 


Figure  13.  PP/P  spectral  ratio  inversion  for  7/6/77 
recorded  at  station  BKS  with  6t*=.67. 

For  the  waveform  comparison,  the  observed 
PP  wave  is  shown  by  the  solid  line,  the 
synthetic  PP  is  shown  by  the  dashed  line. 
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St*  VS  NUMBER  OF  OBSERVATIONS 
P/P P  SPECTRAL  RATIO  METHOD 
PACIFIC  OCEAN 
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Ficjure  14. 


Mistoqram  of  tft*  inversion  results  for  Pacific 
Ocean  bounce  point. 
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ditiun  requires  a  significant  source  anisotropy,  such  as 
takeoff  angles  near  the  nodal  planes. 

The  conclusions  to  be  reached  from  this  early  study  of 

k 

PP/P  inversion  for  6t  is  that  bandwidth  is  the  greatest 

k 

problem.  There  simply  is  not  enough  t  signal  in  the  low 

* 

frequency  crustal  noise  to  permit  clear  estimation  of  t  . 
Adequate  results  will  be  achieved  only  when  short  period  or 
broadband  data  are  brought  to  bear  upon  the  problem.  The 
method  is,  however,  quite  effective,  and  it  will  not  be 

k 

limited  by  assumptions  of  frequency  independent  t  when  high 
frequencies  are  added  to  the  data  set. 

3.2  Receiver  Function  Spectral  Ratios 

Two  data  sets  will  be  used  in  this  section,  rep¬ 
resenting  studies  of  receiver  functions  for  stations  about 
the  Eastern  Kazakh  Test  Site  and  the  Nevada  Test  Site. 
Table  2  lists  the  stations  used  and  their  distance  and 
azimuth  relation  to  the  source  regions. 

An  example  of  an  inversion  is  shown  in  Figure  15  for 
the  ratio  of  CHTO  to  Alice  Springs,  Australia  (ASA)  from  an 
Eastern  Kazakh  source.  In  this  case  the  reference  path  had 

k 

parameters  t  =  .63  sec  and  =  -15  sec;  while  the  path  to 

k  . 

CHTO  has  a  nearly  frequency  independent  t  =  .21  sec.  This 
particular  spectral  ratio  has  a  great  deal  of  character,  and 
very  similar  results  are  found  for  a  wide  range  of  starting 
values  and  weights.  But  examination  of  the  tradeoff  matrix 
shows  that  only  differences  between  the  t  or  x ^  values  are 
tightly  constrained. 
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Figure  15.  t*  inversion  on  receiver  function  spectral  ratio  of 
CHTO  to  ASA  for  an  Eastern  Kazakh  source. 
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Table  2.  Station  Set  for  t* 

Inversion 

Station 

A 

Azimuth 

eastern  Kazakh  Test  Site 

Alice  Springs  (ASA) 

88° 

131° 

ANMO 

95° 

4° 

CHTO 

35° 

145° 

KSKS 

37° 

90° 

PINEDALE 

83° 

J  V 

340° 

Nevada  Test  Site 

BOCO 

51° 

119° 

GRFO 

82° 

31° 

KONO 

74° 

25° 

KSRS 

85° 

314° 

MAJO 

79° 

308° 

NAO 

73° 

24° 

ZOBO 

70° 

130° 
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In  order  to  reduce  the  tradeoff  problems,  several 
spectral  ratios  may  be  put  into  one  inversion.  Using  the 
same  reference  station  for  each  ratio,  an  inversion  on  four 
ratios  can  be  solved  for  the  properties  of  five  paths.  In 
such  a  case  the  tradeoff  between  individual  paths  is  min- 

k 

imal,  leaving  tradeoff  between  t  and  X2  for  each  path.  In 
practice,  if  the  spectral  ratios  have  a  great  deal  t  f  char¬ 
acter,  few  combinations  exist  which  will  satisfy  a  four 
ratio  inversion. 

Figures  16  and  17  show  example  of  four  ratio  inversions 
for  the  two  source  regions  using  reference  stations  ASA  and 
BOCO,  respectively.  Note  that  the  values  found  for  the  CHTO 
ratio  compare  very  well  to  those  shown  in  Figure  15.  Com¬ 
parison  may  also  be  made  between  the  two  values  for  KSRS, 
the  only  station  common  to  the  two  arrays.  If  the  results 
are  taken  at  face  value,  then  the  attenuation  to  KSRS  from 
Nevada  is  slightly  greater  than  from  central  Asia.  In  a 

difficult  comparison,  however,  although  the  paths  to  KONO 

.  * 

and  NAO  are  very  similar,  the  combination  of  St  and  ST2 
show  significantly  different  path  attenuation  from  NTS. 

The  conclusion  must  be  that  the  method  shows  promise, 
but  that  not  all  problems  have  been  solved.  Two  separate 
classes  of  problem  may  be  defined:  inadequate  or  bandlimited 
data,  and  the  intricacies  of  the  inverse  problem.  The  first 
problem  must  be  solved  by  augmenting  the  current  data  sets, 
and  the  second  problem  may  be  approached  by  testing  the 
procedure  under  a  variety  of  conditions.  A  few  such  tests 
will  now  be  discussed. 
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In  order  to  equalize  the  importance  of  each  model 
parameter,  the  weights  should  be  adjusted  such  that  the 
eigenvalues  of  equations  14-16  are  of  equal  size.  Because 

•k 

t  is  on  the  order  of  1.0  sec  and  is  on  the  order  of  0.1 
m  2 

sec,  weights  were  set  at  0.1  and  1.0,  respectively,  where 

layer  weights  reduce  the  freedom  of  the  variable.  Of 

course,  the  weights  may  be  fine-tuned  once  appropriate 

starting  values  have  been  determined. 

With  the  weights  set,  the  behavior  as  a  function  of 

starting  model  may  be  tested.  Table  3  shows  the  result  of 

* 

trying  four  different  starting  models  for  which  the  t  and 

r2  trade  off  to  give  roughly  equivalent  path  attenuation. 

*  ... 
The  tradeoff  between  t  and  i2  at  each  station  is  evident  in 

the  table,  but  the  6t  and  St2  are  consistent  independent  of 

starting  values.  The  path  from  Kazakh  to  Thailand  (CHTO) 

* 

consistently  gives  t  on  the  order  of  .25  sec  and  very 
small  t2.  In  the  only  conflicting  result,  the  squared  error 

•k  . 

increases  by  35%.  The  6 1  between  the  reference  station, 
ASA,  and  ANMO  is  consistently  greater  than  .85,  indicating  a 
path  with  high  attenuation.  Of  the  other  paths,  those  to 
KSRS  and  PINEDALE  show  slightly  higher  attenuation  than  that 
to  ASA,  independent  of  starting  model. 

Table  4  shows  a  test  based  upon  stations  in  the  array 
for  the  NTS  receiver  function  spectral  ratios.  In  these 
tests,  the  constraints  on  the  inverse  are  changed  by  in¬ 
cluding  different  data  in  each  run;  so  the  consistency 
demonstrated  in  the  table  is  very  important.  The  most 


j 

i 
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■k  t 

TABLE  3.  t  Inversion  on  Receiver  Function  Spectral  Ratios 
Variation  with  Parameter  Starting  Values 


Starting  ALICE 

Value  SPRINGS  ANMO  CHTO  KSRS  PINEDALE  ERROR 

_ t*  x2  t*  T2  T  2  T  2  X2 

t*  -  .45  .32  .06  1.17  .07  .22  .00  .36  .05  .39  .05  29.1 

t2  =  .05 

t*  =  .60  .46  .09  1.41  .09  .24  .01  .52  .08  .55  .07  28.8 

t2  =  .08 

t*  =  .80  .68  .12  1.64  .11  .28  .02  .75  .12  .76  .10  30.8 

x2  =  .12 

t*  =  1.00  .83  .17  1.77  .13  .69  .09  .91  .16  .95  .13  42.7 

T  2  =  -15 


SGI -R-81-052 


I* 


42 


* 

TABLE  4.  t  Inversion  on  Receiver  Function  Spectral  Ratios 


BOCO 

Variation  with  Stations  in  Array 

GRFO  KONO  KSRS  MAJO 

NAO 

ERROR 

* 

t 

X2 

t* 

t2 

* 

t 

X2 

* 

t 

X2 

* 

t 

X  2 

* 

t  l 

2 

.91 

.11 

.80  .06 

.57 

.02 

.64 

.04 

.57 

.04 

9.2 

.94 

.10 

.88  .06 

.51 

.01 

.52 

.03 

.69 

.04 

8.2 

1.08 

.10 

1.08  .07 

.56 

.02 

.38 

.00 

.63 

.03 

6.9 

1.16 

.09 

.57 

.00 

.71 

.03 

.55 

.02 

.78 

.03 

8.0 

.91 

.10 

.86  .06 

.49 

.01 

.61 

.03 

.68 

.04 

8.9 
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significant  change  is  caused  by  removing  station  KONO,  but 

examination  of  Figure  17  shows  that  KONO  ratio  has  more 

character  than  the  others.  Thus  it  provides  the  strongest 

constraint  on  the  inversion. 

Changing  the  reference  station  changes  only  one  path  in 

the  inverse,  but  it  changes  the  noise  properties  of  all  of 

the  ratios.  This  in  turn  changes  the  tradeoff  properties  of 

the  inversion,  as  is  demonstrated  by  Figure  18.  In  spite  of 

the  radical  change  in  reference  path  properties,  however, 

the  same  relative  secondary  path  properties  are  still  found 

* 

by  the  inversion.  Table  5  shows  the  fit  between  the  sec¬ 
ondary  stations  for  both  inversion  runs,  and  the  patterns 
are  nearly  exact. 

The  final  test  is  of  the  model  definition  itself. 

Examination  of  equation  (3)  shows  that,  in  the  limit  as  t2 

goes  to  zero,  the  inverse  tangent  reaches  its  maximum  and 
£  ★ 

t  (uj)  =  t  .  This  is  the  commonly  used  constant  Q  assump- 

* 

tion.  Another  common  assumption  is  that  t  =  1.0  sec  at 

low  frequencies  (Lundquist,  et  al.,  1980).  In  that  case, 

*  *  . 

the  frequency  dependence  of  t  is  defined  entirely  by 
Both  of  these  assumptions  reduce  the  number  of  model  para¬ 
meters,  thus  reducing  the  ability  of  the  inversion  to  fit 

the  data.  If  the  quality  of  the  fit  is  significantly  worse, 

* 

then  the  assumption  of  frequency  dependent  t  is  justified. 

Figures  19-22  show  the  parameter  reduction  tests  for 

the  two  arrays,  and  comparison  should  be  made  to  Figures  16 

★ 

and  17  which  show  equivalent  runs  for  t  (f).  The  results 
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TABLE  5.  Inversion  Results  and  6t 
with  Changing  Reference  Station 


BOCO  Reference  ZOBO  Reference 


* 

t 

T  2 

* 

fit 

t 

T2 

* 

fit 

GRFO 

.86 

.06 

.37 

1.09 

.  10 

.32 

KONO 

.49 

.01 

-.12 

.77 

.06 

-.11 

KSRS 

.61 

.03 

-.07 

.88 

.08 

-.06 

NAO 

.68 

.  04 

-.18 

.94 

.08 

-.15 

GRFO 

.86 

.06 

1.09 

.10 

I 


Figure  18.  t*  inversion  on  four  receiver  function  spectral  ratios 


using  reference  station  ZOBO  for  NTS  sources.  Compare 
these  results  with  the  results  from  the  previous  figure. 
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are  summarized  in  Table  6  which  shows  t  ,  t 2  and  the  squared 
error.  The  most  apparent  feature  of  Table  6  is  that  the 
different  assumptions  do  not  produce  compatible  results  in 

common  parameters.  CHTO  and  KONO  have  significantly  smaller 

* 

t  than  the  other  stations  (paths)  in  their  respective 
*  *  *  * 

arrays  when  t  =  t  (f).  But  when  t  =  t  ,  each  has  a  large 
* 

t  relative  to  the  other  stations.  That  is,  conclusions 
based  upon  a  frequency  independent  Q  are  in  direct  conflict 
with  results  based  upon  frequency  dependent  Q.  Similar 

: k 

statements  may  be  made  about  i2  under  the  assumption  t  m  = 

1.0. 

The  question  still  remains  as  to  which  assumption  is 

* 

more  correct.  Based  upon  the  squared  error,  the  t  (f) 
assumption  is  best,  but  that  model  has  extra  parameters  with 
which  to  fit  the  data.  The  difference  in  errors  for  the  NTS 
array  is  not  significant,  but  the  difference  for  the  Eastern 
Kazakh  array  definitely  is.  Once  again,  paths  with  unusual 
attenuation  provide  character  in  the  spectral  ratio  which 
can  resolve  the  choice  of  models.  Comparison  of  Figures  16 
and  19  shows  that  the  constant  t  assumption  fails  to  fit 
the  ratios  at  CHTO  and  ANMO .  Within  the  resolution  of  the 
method,  the  constant  t  assumption  is  not  an  adequate  model 
of  path  attenuation. 

The  definitive  answer  must  wait,  however,  for  accurate 
instrument  gains  and/or  broadband  spectral  ratios.  The 
three  tests  for  ANMO  (Figures  16,  19,  21),  for  instance, 
differ  greatly  at  the  low  frequency  end,  so  the  addition  of 
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TABLE 

6.  Inversion  Results 

and  Squared 

Error 

Variation 

with  Model 

Definition 

* 

k 

* 

* 

* 

t  =t 

(f) 

t  =t 

t  = 

1.0 

m 

m 

* 

* 

* 

t 

1  2 

t 

1 2 

t 

X  2 

ASA 

.61 

.12 

.59 

.  001 

1.0 

.05 

ANMO 

1.67 

.11 

.80 

.001 

1.0 

.03 

CHTO 

.20 

.00 

.69 

.001 

1.0 

.03 

KSRS 

.55 

.09 

.63 

.001 

1.0 

.04 

Pinedale 

.70 

.09 

.65 

.001 

1.0 

.04 

ERROR 

28 

.0 

41. 

4 

32 

.4 

BOCO 

.91 

.10 

.44 

.001 

1.0 

.04 

GRFO 

.86 

.06 

.60 

.001 

1.0 

.02 

KONO 

.49 

.01 

.70 

.001 

1.0 

.03 

KSRS 

.61 

.03 

.63 

.001 

1.0 

.02 

NAO 

.68 

.04 

.63 

.001 

1.0 

.02 

ERROR 

8. 

9 

10. 

5 

9. 

5 
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long  period  contraints  can  resolve  the  choice  of  models.  If 
the  relative  gains  shown  here  are  correct,  then  the  ratio 
should  go  through  unity  at  low  frequency.  In  that  case,  a 
curve  with  two  inflections  is  required  to  fit  the  ANMO 
spectral  ratios,  and  neither  of  the  other  assumptions  is 
adequate . 


I 


I* 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 

Both  PP/P  and  receiver  function  spectral  ratios  have 

: 

been  interpreted  in  terms  of  t  (f),  and  both  tests  have 
found  that  inadequate  data  sets  limit  the  resolution  of 
model  parameters.  In  each  case  significant  noise  within  the 
passband  of  the  ratios  was  due  to  crustal  transfer  func¬ 
tions.  Among  other  conclusions  to  be  made  here  is  the  fact 
that  crustal  effects  are  very  significant,  much  more  im¬ 
portant  than  might  be  estimated  via  a  Haskell-Thompson 
formulation . 

Tradeoff  problems  also  limit  the  accuracy  of  the  in- 
version  technique.  t  and  t2  may  trade  off  to  give  roughly 
equivalent  apparent  attenuation  over  any  narrow  bandwidth, 
and  the  parameters  for  two  paths  can  trade  off  in  a  spectral 

ratio.  Thus  differences  are  more  correctly  estimated  than 

* 

absolute  values.  The  PP/P  ratios  determine  a  6t  which  can 
be  associated  with  the  bounce  point,  and  the  most  reliable 

: k 

result  for  central  Asia  gives  t  =0.8  sec  for  the  upper 
mantle  beneath  the  Eastern  Kazakh  test  site. 

Receiver  function  spectral  ratios  can  determine  rel¬ 
ative  path  properties  for  paths  from  a  common  source  region. 
Considering  paths  about  the  Eastern  Kazakh  test  site, 
station  CHTO  sees  very  little  attenuation,  and  ANMO  sees  a 
very  high  attenuation.  The  6t  between  these  two  stations 

is  more  than  1.2  sec.  Stations  KSRS,  Pinedale  and  Alice 

* 

Springs  are  intermediate,  averaging  6t  =.4  greater  than 
CHTO.  For  paths  about  NTS,  ZOBO  sees  a  very  low  path 
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attenuation,  and  KONO  has  somewhat  lower  attenuation  than 
KSRS  or  NAO.  The  6t  -.2  sec  between  KONO  and  NAO  is  of 
interest  because  of  the  proximity  of  the  two  stations. 

k 

It  is  of  interest  to  note  that  t  values  found  by  this 

inversion  technique  cluster  in  the  range  .4  to  1.0  sec,  and 

t2  is  never  found  to  be  greater  than  .18  sec.  Though  this 

is  somewhat  dependent  upon  starting  models,  it  is  clear  that 
* 

t  values  do  not  tend  toward  values  greater  than  1.0. 

Perhaps  the  most  important  conclusion  sought  by  this 
research  is  a  definitive  statement  as  to  the  frequency 
dependence  of  Q.  Within  the  limitations  of  the  data  sets 
and  the  method,  it  appears  that  the  assumption  of  constant  Q 

k 

(or  constant  t  )  does  not  provide  an  adequate  fit  to  the 
data  sets.  An  unequivocal  conclusion,  however,  will  require 
more  work,  data  sets  with  greater  bandwidth,  and  accurate 
instrument  calibrations. 
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